Evolutionary ecology of host competence after a chytrid outbreak in a naive amphibian community.

Naive multi-host communities include species that may differentially maintain, transmit and amplify novel pathogens; therefore, we expect species to fill distinct roles during infectious disease emergence. Characterizing these roles in wildlife communities is challenging because most disease emergence events are unpredictable. Here, we used field-collected data to investigate how species-specific attributes influenced the degree of exposure, probability of infection, and pathogen intensity, during the emergence of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in a highly diverse tropical amphibian community. Our findings confirmed that ecological traits commonly evaluated as correlates of decline were positively associated with infection prevalence and intensity at the species level during the outbreak. We identified key hosts that disproportionally contributed to transmission dynamics in this community and found a signature of phylogenetic history in disease responses associated with increased pathogen exposure via shared life-history traits. Our findings establish a framework that could be applied in conservation efforts to identify key species driving disease dynamics under enzootics before reintroducing amphibians back into their original communities. Reintroductions of supersensitive hosts that are unable to overcome infections will limit the success of conservation programmes by amplifying the disease at the community level. This article is part of the theme issue 'Amphibian immunity: stress, disease and ecoimmunology'.

Host immune responses to enzootic and invasive pathogen lineages vary in magnitude, timing, and efficacy.

Infectious diseases of wildlife continue to pose a threat to biodiversity worldwide, yet pathogens are far from uniform in virulence or host disease outcome. Within the same pathogen species, virulence can vary considerably depending on strain or lineage, in turn eliciting variable host responses. One pathogen that has caused extensive biodiversity loss is the amphibian-killing fungus, Batrachochytrium dendrobatidis (Bd), which is comprised of a globally widespread hypervirulent lineage (Bd-GPL), and multiple geographically restricted, enzootic lineages. Whereas host immunogenomic responses to Bd-GPL have been characterized in a number of amphibian species, immunogenomic responses to geographically restricted, enzootic Bd lineages are less clear. To examine lineage-specific host immune responses to Bd, we exposed a species of pumpkin toadlet, Brachycephalus pitanga, which is endemic to Brazil's Southern Atlantic Forest, to either the Bd-GPL or the enzootic Bd-Asia-2/Brazil (hereafter Bd-Brazil) lineage. Using temporal samples from early, mid, and late infection stages, we quantified functional immunogenomic responses over the course of infection using differential gene expression tests and coexpression network analyses. Host immune responses varied significantly with Bd lineage. Relative to controls, toadlet responses to Bd-Brazil were weak at early infection (25 genes significantly differentially expressed), peaked by mid-stage infection (414 genes), and were nearly fully resolved by late-stage infection (nine genes). In contrast, responses to Bd-GPL were magnified and delayed; toadlets significantly differentially expressed 111 genes early, 87 genes at mid-stage infection, and 726 genes by late-stage infection relative to controls. Given that infection intensity did not vary between mid- and late-stage disease in either Bd-Brazil or Bd-GPL treatments, this suggests that pumpkin toadlets may be at least partially tolerant to the enzootic Bd-Brazil lineage. In contrast, late-stage immune activation against Bd-GPL was consistent with immune dysregulation previously observed in other species. Our results demonstrate that both the timing of immune response and the particular immune pathways activated are specific to Bd lineage. Within regions where multiple Bd lineages co-occur, and given continued global Bd movement, these differential host responses may influence not only individual disease outcome, but transmission dynamics at the population and community levels.

Environment rather than character displacement explains call evolution in glassfrogs.

The acoustic adaptation hypothesis (AAH) and ecological character displacement (ECD) are two potential mechanisms shaping call evolution that can predict opposite trends for the differentiation of signals. Under AAH, signals evolve to minimize environmental degradation and maximize detection against background noise, predicting call homogenization in similar habitats due to environmental constraints on signals. In contrast, ECD predicts greater differences in call traits of closely related taxa in sympatry because of selection against acoustic interference. We used comparative phylogenetic analyses to test the strength of these two selective mechanisms on the evolution of advertisement calls in glassfrogs, a highly diverse family of neotropical anurans. We found that, overall, acoustic adaptation to the environment may outweigh effects of species interactions. As expected under the AAH, temporal call parameters are correlated with vegetation density, but spectral call parameters had an unexpected inverse correlation with vegetation density, as well as an unexpected correlation with temperature. We detected call convergence among co-occurring species and also across multiple populations from the same species in different glassfrogs communities. Our results indicate that call convergence is common in glassfrogs, likely due to habitat filtering, while character displacement is relatively rare, suggesting that costs of signal similarity among related species may not drive divergent selection in all systems.

Hybridization despite elaborate courtship behavior and female choice in Neotropical tree frogs.

The mechanisms of hybridization can be elucidated by analyzing genotypes as well as phenotypes that could act as premating barriers, as the reproductive interactions among heterospecifics can alter the evolutionary history of species. In frogs, hybrids typically occur among species that reproduce explosively (in dense aggregations) with few opportunities for mate selection but are rare in species with elaborate courtship behaviors that may prevent erroneous mating. Using 21 microsatellite markers, we examined hybridization in the prolonged-breeding tree frogs Bokermannohyla ibitiguara and B. sazimai sampled within a contact zone in the Brazilian savanna (72 tadpoles; 74 adults). We also compared acoustic and morphological data. We confirmed both parental species genetically; STRUCTURE results confirmed 14 hybrids, 11 of which were second-generation according to NEWHYBRIDS, all with intermediate values of genetic dissimilarities compared to the parentals. Morphological and acoustic analyses revealed that hybrids showed variable but not necessarily intermediate phenotypes. Moreover, 2 hybrids exhibited call types different from parentals. The reproduction of B. ibitiguara involves territorial and aggressive males, elaborate courtships with acoustic and tactile stimuli, choosy females, and opportunistic strategies. Our study uncovers a rare case of viable hybridization among closely related frogs with such a combination of complex courtship behaviors and mate choice. We discuss the likely directionality and mechanisms behind this phenomenon, and highlight the importance of investigating hybridization even in species that show elaborate reproduction and female choice to advance our understanding of animal diversification.

Cryptic lineages and standing genetic variation across independent cane toad introductions.

Widespread introduced species can be leveraged to investigate the genetic, ecological and adaptive processes underlying rapid evolution and range expansion, particularly the contributions of genetic diversity to adaptation. Rhinella marina, the cane toad, has been a focus of invasion biology for decades in Australia. However, their introduction history in North America is less clear. Here, we investigated the roles of introduction history and genetic diversity in establishment success of cane toads across their introduced range. We used reduced representation sequencing (ddRAD) to obtain 34,000 SNPs from 247 toads in native (French Guiana, Guyana, Ecuador, Panama, Texas) and introduced (Bermuda, southern Florida, northern Florida, Hawai'i, Puerto Rico) populations. Unlike all other cane toad introductions, we found that Florida populations were more closely related to native Central American lineages (R. horribilis), than to native Southern American lineages (R. marina). Furthermore, we found high levels of diversity and population structure in the native range, corroborating suggestions that R. marina is a species complex. We also found that introduced populations exhibit only slightly lower genetic diversity than native populations. Together with demographic analyses, this indicates founding populations of toads in Florida were larger than previously reported. Lastly, within R. marina, only one of 245 putatively adaptive SNPs showed fixed differences between native and introduced ranges, suggesting that putative selection in these introduced populations is based upon existing genetic variation. Our findings highlight the importance of genetic sequencing in understanding biological introductions and hint at the role of standing genetic variation in range expansion.

Temporal and spatial diversification along the Amazonia-Cerrado transition in Neotropical treefrogs of the Boana albopunctata species group.

Despite extensive research on biodiversity in Neotropical forests, biodiversity in seasonally dry, open biomes in South America has been underestimated until recently. We leverage a widespread group, Boana albopunctata, to uncover cryptic lineages and investigate the timing of diversification in Neotropical anurans with a focus on dry diagonal biomes (Cerrado, Caatinga and Chaco) and the ecotone between Amazonia and the Cerrado. We inferred a multilocus phylogeny of the B. albopunctata species group that includes 15 of 18 described species, recovered two cryptic species, and reconstructed the timing of diversification among species distributed across multiple South American biomes. One new potential species (B. aff. steinbachi), sampled in the Amazonian state of Acre, clustered within the B. calcara-fasciata species complex and is close to B. steinbachi. A second putative new species (B. aff. multifasciata), sampled in the Amazonia-Cerrado ecotone, is closely related to B. multifasciata. Lastly, we place a recently identified Cerrado lineage (B. aff. albopuncata) into the B. albopunctata species group phylogeny for the first time. Our ancestral range reconstruction showed that species in the B. albopuctata group likely dispersed from Amazonia-Cerrado into the dry-diagonal and Atlantic Forest. Intraspecies demography showed, for both B. raniceps and B. albopunctata, signs of rapid expansion across the dry diagonal. Similarly, for one clade of B. multifasciata, our analyses support an invasion of the Cerrado from Amazonia, followed by a rapid expansion across the open diagonal biomes. Thus, our study recovers several recent divergences along the Amazonia-Cerrado ecotone in northern Brazil. Tectonic uplift and erosion in the late Miocene and climate oscillations in the Pleistocene corresponded with estimated divergence times in the dry diagonal and Amazonia-Cerrado ecotone. Our study highlights the importance of these threatened open formations in the generation of biodiversity in the Neotropics.

Endemic Lineages of Batrachochytrium dendrobatidis Are Associated With Reduced Chytridiomycosis-Induced Mortality in Amphibians: Evidence From a Meta-Analysis of Experimental Infection Studies.

Emerging infectious wildlife diseases have caused devastating declines, particularly when pathogens have been introduced in naïve host populations. The outcome of disease emergence in any host population will be dictated by a series of factors including pathogen virulence, host susceptibility, and prior opportunity for coevolution between hosts and pathogens. Historical coevolution can lead to increased resistance in hosts and/or reduced virulence in endemic pathogens that allows stable persistence of host and pathogen populations. Adaptive coevolution may also occur on relatively short time scales following introduction of a novel pathogen. Here, we performed a meta-analysis of multi-strain Batrachochytrium dendrobatidis (Bd) infection experiments to test whether: (1) amphibian hosts exhibit lower mortality rates when infected with strains belonging to endemic Bd lineages relative to the Global Panzootic Lineage (Bd-GPL), hypothetically owing to long co-evolutionary histories between endemic Bd lineages and their amphibian hosts; and (2) amphibians exhibit lower mortality rates when infected with local Bd-GPL strains compared with non-local Bd-GPL strains, hypothetically owing to recent selection for tolerance or resistance to local Bd-GPL strains. We found that in a majority of cases, amphibians in endemic Bd treatments experienced reduced mortality relative to those in Bd-GPL treatments. Hosts presumed to have historically coexisted with endemic Bd did not show reduced mortality to Bd-GPL compared with hosts that have not historically coexisted with endemic Bd. Finally, we detected no overall difference in amphibian mortality between local and non-local Bd-GPL treatments. Taken together, our results suggest that long-term historical coexistence is associated with less disease-induced mortality potentially due to hypovirulence in endemic Bd lineages, and that more recent coexistence between amphibians and Bd-GPL has not yet resulted in reduced host susceptibility or pathogen virulence. This corroborates previous findings that Bd-GPL introduced via the global amphibian trade has a high capacity for causing disease-induced mortality.

Meta-analysis of Gender Performance Gaps in Undergraduate Natural Science Courses.

To investigate patterns of gender-based performance gaps, we conducted a meta-analysis of published studies and unpublished data collected across 169 undergraduate biology and chemistry courses. While we did not detect an overall gender gap in performance, heterogeneity analyses suggested further analysis was warranted, so we investigated whether attributes of the learning environment impacted performance disparities on the basis of gender. Several factors moderated performance differences, including class size, assessment type, and pedagogy. Specifically, we found evidence that larger classes, reliance on exams, and undisrupted, traditional lecture were associated with lower grades for women. We discuss our results in the context of natural science courses and conclude by making recommendations for instructional practices and future research to promote gender equity.

Implementing team-based learning in the life sciences: A case study in an online introductory level evolution and biodiversity course.

Team-Based Learning (TBL) is a pedagogical tool that has great potential to develop student engagement, accountability, and equity in the online classroom. TBL is rooted in evidence-based educational theories and practices that underlie many active learning approaches such as self-testing, team discussion, and application of knowledge. The use of these approaches is associated with better student performance, retention, and sense of belonging in the classroom, aspects that are often reported to be especially lacking in online courses. Here, we describe how we implemented TBL in a face-to-face and an online introductory level evolution and biodiversity course. We implemented TBL in the face-to-face course (~200 students) starting in 2018 and in the online course (~30 students) starting in the summer of 2019. We used several online applications to facilitate the transition to an online platform such as Simbio, Slack, VoiceThread, Articulate 360, and Teammates. Our experiences using TBL approaches in the online course have been rewarding, and students are engaged and accountable for their learning and performed well in the course. Our goal is to provide an example of how we designed a life science course using TBL approaches and transitioned the course to an online environment. With the current switch to remote instruction and online learning, we recommend the use of TBL as a course design approach that can improve the students' online learning experience.

Temperature dependence of metabolic rate in tropical and temperate aquatic insects: Support for the Climate Variability Hypothesis in mayflies but not stoneflies.

A fundamental gap in climate change vulnerability research is an understanding of the relative thermal sensitivity of ectotherms. Aquatic insects are vital to stream ecosystem function and biodiversity but insufficiently studied with respect to their thermal physiology. With global temperatures rising at an unprecedented rate, it is imperative that we know how aquatic insects respond to increasing temperature and whether these responses vary among taxa, latitudes, and elevations. We evaluated the thermal sensitivity of standard metabolic rate in stream-dwelling baetid mayflies and perlid stoneflies across a ~2,000 m elevation gradient in the temperate Rocky Mountains in Colorado, USA, and the tropical Andes in Napo, Ecuador. We used temperature-controlled water baths and microrespirometry to estimate changes in oxygen consumption. Tropical mayflies generally exhibited greater thermal sensitivity in metabolism compared to temperate mayflies; tropical mayfly metabolic rates increased more rapidly with temperature and the insects more frequently exhibited behavioral signs of thermal stress. By contrast, temperate and tropical stoneflies did not clearly differ. Varied responses to temperature among baetid mayflies and perlid stoneflies may reflect differences in evolutionary history or ecological roles as herbivores and predators, respectively. Our results show that there is physiological variation across elevations and species and that low-elevation tropical mayflies may be especially imperiled by climate warming. Given such variation among species, broad generalizations about the vulnerability of tropical ectotherms should be made more cautiously.

Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil's Atlantic Forest treefrogs.

In Brazil's Atlantic Forest (AF) biodiversity conservation is of key importance since the fungal pathogen Batrachochytrium dendrobatidis (Bd) has led to the rapid loss of amphibian populations here and worldwide. The impact of Bd on amphibians is determined by the host's immune system, of which the skin microbiome is a critical component. The richness and diversity of such cutaneous bacterial communities are known to be shaped by abiotic factors which thus may indirectly modulate host susceptibility to Bd. This study aimed to contribute to understanding the environment-host-pathogen interaction determining skin bacterial communities in 819 treefrogs (Anura: Hylidae and Phyllomedusidae) from 71 species sampled across the AF. We investigated whether abiotic factors influence the bacterial community richness and structure on the amphibian skin. We further tested for an association between skin bacterial community structure and Bd co-occurrence. Our data revealed that temperature, precipitation, and elevation consistently correlate with richness and diversity of the skin microbiome and also predict Bd infection status. Surprisingly, our data suggest a weak but significant positive correlation of Bd infection intensity and bacterial richness. We highlight the prospect of future experimental studies on the impact of changing environmental conditions associated with global change on environment-host-pathogen interactions in the AF.

Unexpected reproductive fidelity in a polygynous frog.

Polygynous mating systems with group fidelity are a common animal organization, typically consisting of multiple females in a mated group with a single male for an extended period (sometimes referred to as harem polygyny). Single-male polygyny with reproductive fidelity occurs in invertebrates, bony fishes, and some tetrapods, such as lizards, mammals, and birds. In amphibians, reproductive fidelity in polygynous groups is not fully demonstrated. Combining data on larval development, molecular paternity assignment, and in situ behavioral observations, we reveal high fidelity during a prolonged breeding season in a Neotropical polygynous frog. Males dominate scarce breeding sites, guarding offspring, and mating exclusively with multiple females that exhibit dominance rank. This system likely evolved in response to intense competition for breeding sites and intrasexual competition for mates.

Incapacitating effects of fungal coinfection in a novel pathogen system.

As globalization lowers geographic barriers to movement, coinfection with novel and enzootic pathogens is increasingly likely. Novel and enzootic pathogens can interact synergistically or antagonistically, leading to increased or decreased disease severity. Here we examine host immune responses to coinfection with two closely related fungal pathogens: Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Both pathogens have had detrimental effects on amphibian populations, with Bd now largely enzootic, while Bsal is currently spreading and causing epizootics. Recent experimental work revealed that newts coinfected with Bd and Bsal had significantly higher mortality than those infected with either pathogen alone. Here we characterize host immunogenomic responses to chytrid coinfection relative to single infection. Across several classes of immune genes including pattern recognition receptors, cytokines, and MHC, coinfected host gene expression was weakly upregulated or comparable to that seen in single Bd infection, but significantly decreased when compared to Bsal infection. Combined with strong complement pathway downregulation and keratin upregulation, these results indicate that coinfection with Bd and Bsal compromises immune responses active against Bsal alone. As Bsal continues to invade naïve habitats where Bd is enzootic, coinfection will be increasingly common. If other Bd-susceptible species in the region have similar responses, interactions between the two pathogens could cause severe population and community-level declines.

Response to Comment on "Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity".

Lambert et al question our retrospective and holistic epidemiological assessment of the role of chytridiomycosis in amphibian declines. Their alternative assessment is narrow and provides an incomplete evaluation of evidence. Adopting this approach limits understanding of infectious disease impacts and hampers conservation efforts. We reaffirm that our study provides unambiguous evidence that chytridiomycosis has affected at least 501 amphibian species.

Male-male competition and repeated evolution of terrestrial breeding in Atlantic Coastal Forest frogs.

Terrestrial breeding is a derived condition in frogs, with multiple transitions from an aquatic ancestor. Shifts in reproductive mode often involve changes in habitat use, and these are typically associated with diversification in body plans, with repeated transitions imposing similar selective pressures. We examine the diversification of reproductive modes, male and female body sizes, and sexual size dimorphism (SSD) in the Neotropical frog genera Cycloramphus and Zachaenus, both endemic to the Atlantic rainforest of Brazil. Species in this clade either breed in rocky streams (saxicolous) or in terrestrial environments, allowing us to investigate reproductive habitat shifts. We constructed a multilocus molecular phylogeny and inferred evolutionary histories of reproductive habitats, body sizes, and SSD. The common ancestor was small, saxicolous, and had low SSD. Terrestrial breeding evolved independently three times and we found a significant association between reproductive habitat and SSD, with shifts to terrestrial breeding evolving in correlation with decreases in male body size, but not female body size. Terrestrial breeding increases the availability of breeding sites and results in concealment of amplexus, egg-laying, and parental care, therefore reducing male-male competition at all stages of reproduction. We conclude that correlated evolution of terrestrial reproduction and small males is due to release from intense male-male competition that is typical of exposed saxicolous breeding.

Thermal cues drive plasticity of desiccation resistance in montane salamanders with implications for climate change.

Organisms rely upon external cues to avoid detrimental conditions during environmental change. Rapid water loss, or desiccation, is a universal threat for terrestrial plants and animals, especially under climate change, but the cues that facilitate plastic responses to avoid desiccation are unclear. We integrate acclimation experiments with gene expression analyses to identify the cues that regulate resistance to water loss at the physiological and regulatory level in a montane salamander (Plethodon metcalfi). Here we show that temperature is an important cue for developing a desiccation-resistant phenotype and might act as a reliable cue for organisms across the globe. Gene expression analyses consistently identify regulation of stem cell differentiation and embryonic development of vasculature. The temperature-sensitive blood vessel development suggests that salamanders regulate water loss through the regression and regeneration of capillary beds in the skin, indicating that tissue regeneration may be used for physiological purposes beyond replacing lost limbs.

Close-kin mating, but not inbred parents, reduces hatching rates and offspring quality in a threatened tortoise.

Inbreeding depression, the reduction in fitness due to mating of related individuals, is of particular conservation concern in species with small, isolated populations. Although inbreeding depression is widespread in natural populations, long-lived species may be buffered from its effects during population declines due to long generation times and thus are less likely to have evolved mechanisms of inbreeding avoidance than species with shorter generation times. However, empirical evidence of the consequences of inbreeding in threatened, long-lived species is limited. In this study, we leverage a well-studied population of gopher tortoises, Gopherus polyphemus, to examine the role of inbreeding depression and the potential for behavioural inbreeding avoidance in a natural population of a long-lived species. We tested the hypothesis that increased parental inbreeding leads to reduced hatching rates and offspring quality. Additionally, we tested for evidence of inbreeding avoidance. We found that high parental relatedness results in offspring with lower quality and that high parental relatedness is correlated with reduced hatching success. However, we found that hatching success and offspring quality increase with maternal inbreeding, likely due to highly inbred females mating with more distantly related males. We did not find evidence for inbreeding avoidance in males and outbred females, suggesting sex-specific evolutionary trade-offs may have driven the evolution of mating behaviour. Our results demonstrate inbreeding depression in a long-lived species and that the evolution of inbreeding avoidance is shaped by multiple selective forces.

Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity.

Anthropogenic trade and development have broken down dispersal barriers, facilitating the spread of diseases that threaten Earth's biodiversity. We present a global, quantitative assessment of the amphibian chytridiomycosis panzootic, one of the most impactful examples of disease spread, and demonstrate its role in the decline of at least 501 amphibian species over the past half-century, including 90 presumed extinctions. The effects of chytridiomycosis have been greatest in large-bodied, range-restricted anurans in wet climates in the Americas and Australia. Declines peaked in the 1980s, and only 12% of declined species show signs of recovery, whereas 39% are experiencing ongoing decline. There is risk of further chytridiomycosis outbreaks in new areas. The chytridiomycosis panzootic represents the greatest recorded loss of biodiversity attributable to a disease.

Community richness of amphibian skin bacteria correlates with bioclimate at the global scale.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analysed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions compared with environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host's macro-environment mediate microbial diversity.